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a. Introduction. — The inner boundary condition imposed by the event horizon requires that accretion 
onto astrophysical black holes may manifest the transonic behavior in general [1—3]. For low angular momentum 
inviscid axisymmetric flow under the influence of the post-Newtonian pseudo-Schwarzschild black hole potentials 
multi-transonicity may be realized in the phase-portrait of the stationary accretion solutions- leading to the possibility 
of the formation of steady standing shock in between the saddle type outer sonic point and the center type middle 
critical point [1, 4-13]. 

The inviscid flow assumption, however, requires that the accretion flow is characterized by reasonably low angular 
momentum. Such angular momentum distribution should necessarily be sub-Keplerian. Also, such flow is expected to 
possess ab-initio non-zero value of the advective velocity along the equatorial plane of the flow. Such a configuration 
may perhaps be better described by a quasi-spherical flow rather than a standard thin disc structure with radius 
dependent flow thickness. A conical wedge shaped geometrical configuration may be the appropriate description 
for the axisymmetric flow profile under consideration and such geometry may be a better alternative of the flow in 
hydrostatic equilibrium in the vertical direction (hereafter flow in vertical equilibrium in brief). Results obtained so 
far in the literature for inviscid axisymmetric black hole accretion in vertical equilibrium are, thus, to be compared 
with the results obtained using a more general theoretical set up capable of studying various features of the low 
angular momentum inviscid axisymmetric accretion in all three different flow geometries usually considered in the 
relevant literature, viz., accretion flow with constant thickness, conical wedge shaped flow, and accretion in vertical 
equilibrium (Ref [14] and references therein). In most of the works dealing with the multi-transonic accretion flow 
deal with the stationary integral flow solutions. Transient phenomena are, however, not quite rare in large scale 
astrophysical set up. Hence, it is important to establish that the information obtained from the stationary solutions 
are reliable by ensuring that the steady states are stable for black hole accretion - at least for astrophysically relevant 
time scales. Such a task may be accomplished by perturbing the set of equations describing the steady state accretion 
flow and by demonstrating that such perturbations do not diverge with time. For the pseudo-Schwarzschild black 
hole accretion, it has recently been argued [14] that the stationary inviscid flow described by a barotropic equation 
of state, are stable under linear perturbation of the mass accretion rate. Such observations remain valid for all three 
geometrical flow configurations described in previous paragraphs. Existence of certain mapping of the stationary 
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transonic astrophysical flows (in strong gravity) with certain aspects of the analogue gravity (see, e.g., [15-17] for a 
comprehensive discussion about the analogue gravity phenomena in general) has recently been established [18] for 
pseudo-Schwarzschild axisymmetric accretion onto a non-rotating black hole for different geometric profiles of the 
accreting matter. The corresponding acoustic geometry was realized by studying the transonic properties of the 
stationary accretion solutions, and the numerical value of the associated acoustic surface gravity k was estimated. Off 
late, the conjecture that most (if not all) of the astrophysical black holes possess non-zero values of the spin angular 
momentum (a ^ 0, where a is the Kerr parameter) has gained widespread acceptance [19-39]. It is thus essential 
to understand how the black hole spin influences the salient features of the multi-transonic stationary integral flow 
solutions as well as to check whether the linear perturbation scheme gets altered with the variation of a. Which, in 
turn, would highlight how the space-time structure of the associated acoustic geometry depends on the rotation of 
the black hole. Complete general relativistic treatment to address such issues is indeed non-trivial. Although the 
stationary multi-transonic solutions are possible to obtain for flows in the Kerr metric [40] , full space-time dependent 
stability analysis of the transonic solutions does not seem to be amenable within the framework of the general 
relativity, especially when attempts are made to provide a complete analytical solution. 

Motivated by the aforementioned arguments, in the present work we study low angular momentum inviscid 
axisymmetric hydrodynamic accretion onto a rotating black hole using post-Newtonian pseudo-Kerr potentials. 
Several such potentials are available in the literature [41-49]. Among all potentials proposed, we however, pick up 
the one introduced by Artemova et. al [42] since it has the algebraically simplest form. The primary aim behind 
the introduction of a pseudo-Schwarzschild / Kerr black hole potential is to compromise between the convenience 
of handing of a Newtonian description of gravity and the space time structure described by complicated general 
relativistic calculations. Introduction of such potentials should allow one to investigate the basic physical process 
which occur during the accretion phenomena within the framework of a Newtonian construct, hereby avoiding 
the formal general relativistic description of the dynamics of the inflowing matter. The intrinsic simplicity of the 
algebraic presentation of a pseudo-potential is of paramount importance to let that potential qualify as the useful 
one. No particular pseudo-potential stands alone as a significantly superior one (compared to the rest of such 
potentials) to mimic all the properties of the full general relativistic Kerr metric - especially while describing the 
dynamics of the multi-transonic flow - the potential with the simplest algebraic form will be of significant interest for 
the aforementioned purpose. 

The viscous transport of angular momentum is not taken into account in this paper. Forty one years after the 
introduction of the Standard Disc model [50, 51], exact modelling of viscous transonic accretion flow incorporating 
proper heating and cooling effects is still an arduous task. Also, as will be demonstrated later in this work, 
construction of acoustic geometry becomes practically impossible for dissipative flow since the viscous effect may 
destroy the Lorenz in variance. Large radial (advective) velocity close to the horizon ensures that the viscous time is 
larger than the dynamical infall time scale. Large radial velocity even at the larger distance is due to the fact that 
the rotational energy of the axisymmetric flow remains relatively low [52-54], Aforementioned arguments indicate 
that the assumption of inviscid flow is not unjustified for low angular momentum flow with initial non-vanishing 
advective velocity profile. 

For large scale relativistic astrophysical flows, however, the effect of the viscosity and magnetic field may not always 
be neglected. As a consequence, dissipative mechanism (through comptonization, bremsstrahlung or synchrotron 
processes) may become significant to influence the overall flow dynamics and hence the flow may be subjected to the 
turbulent instabilities. Non-linear perturbation may also be of considerable importance to determine the stability 
criteria of the flow. Velikov-Cliandrasekhar [55] and / or Balbus - Hawley [56] type magneto-rotational instabilities 
may be observed close to the black hole event horizon. These additional complexities clearly bring the system to the 
risk of complete destruction of the Lorenz symmetries and hence the analytical modelling of the flow profile as well 
as of the sonic geometry embedded into it becomes impossible - one has to take recourse to the large scale numerical 
simulation to handle such flow profile, which, clearly, is beyond the scope of this work. Also to mention in this aspect 
that the axisymmetric configuration of the flow is assumed a-priori and no misaligned flow structure may be handled 
using the formalism presented in this paper. 

In what follows, we write down the Euler and the continuity equation for matter inflow considered in this 
work. Both polytropic and isothermal flow will be analyzed for three different geometric configurations of the flow. 
These two equations are then integrated to obtain the corresponding first integrals of motion for the respective 
flow configurations. The critical point conditions are derived and the parameter space (spanned by various initial 
boundary conditions governing the flow) for the multi-transonic accretion are identified. The regions of the parameter 
space for which a steady standing stock may form are determined. A linear eigenvalue scheme to identify the nature 
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(whether they are of saddle or of center type) of the critical points is developed. Such a scheme helps to understand 
the nature of the two dimensional Mach number vs. radial distance phase portrait from an algebraic point of view - 
without explicitly evaluating the integral solutions of the governing flow equations. Formation of stationary shock is 
discussed and the parameter sub-space corresponding to the shocked flow is identified. Such a parameter sub-space 
for various flow configurations are compared to understand how the geometric configuration of accreting matter 
influences the shock formation mechanism. The spin dependence of the salient features of the multi-transonic shocked 
flow is investigated. 

Once the properties of the stationary solutions are fully understood, a time dependent linear perturbation analysis 
of the stationary integral flow solutions is performed for a generic form of the pseudo-Kerr black hole potential 
to check whether such perturbations remain non-divergent, to estimate the stability of the steady state accretion 
solutions. It is further demonstrated that a space time dependent acoustic metric can be constructed to describe the 
propagation of the aforementioned linear perturbation and the corresponding pseudo-Riemannian sonic manifold is 
studied. Finally, the black hole spin dependence of the corresponding acoustic surface gravity is studied in somewhat 
detail. 

Hereafter, the radial distance (measured along the equatorial plane of the flow) and any characteristic velocity will 


be scaled in units of 


2 GMbh 


and by c, respectively. Where G, Mbh and c are the universal gravitational constant, 


mass of the black hole considered, and the velocity of light in vacuum, respectively. We adopt natural unit by setting 
G = c = 1. 

The expression for the free fall acceleration as provided by Artemova et-al [42] is 


/ = - 


r 2 ~P(r — r+)^ 


( 1 ) 


The potential corresponding to the above acceleration diverges at the event horizon. 
In the above expression 


r+ = 1 + y/l + a 2 

(2) 

p = Din _ i 

(3) 

r+ 


Tin = 3 + Z 2 — [(3 — ^l)(3 + Zi + 2 .Z 2)] 1 / 2 

(4) 


Z x = 1 + (1 - a 2 ) 1/3 [(l + a) 1 / 3 + (1 - a) 1/3 ] 


(5) 


Z 2 = (3a 2 + Z 2 ) 1 / 2 

Explicit expression for the associated black hole potentials is given as 


r P-2 

r 2 

+ 


r + -/3r + 

r+~pr + 


(j — r+)P 2 


( 6 ) 


( 7 ) 


The Euler and the continuity equation for a generalized pseudo-Kerr potential <f> can be expressed as 
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du 1 dP dd>(r ) 

U d^ + -p^ + — = ° 


( 8 ) 


d[u P rh{r)\ = Q (g) 

dr 

In the present paper, the polytropic as well as the isothermal accretion will be considered. The equation of state 
corresponding to a polytropic flow is p = Kp 1 , p , p and 7 being the pressure, mass, density and the adiabatic index 
( 7 = C p /C v , C p and C v are the specific heats at constant pressure and volume, respectively), respectively. The 
proportionality constant K may be interpreted as a measure of the entropy content of the flow since entropy per 
particle of an ensemble may be expressed as [57], 

1 'Y 

a = -log A' H-- + constant (10) 

7 — 1 7 — 1 


Where the value of the constant depends on the chemical composition of the matter. 


For the isothermal flow, the equation of state is expressed as 


P = 


pRT 

P 


= C 2 sP 


( 11 ) 


Where C s is the position independent (since C s oc T 1 / 2 ) isothermal sound speed. R and p, are the universal gas 
constant and the mean molecular weight, respectively. T is the proton temperature. 


One needs to solve the time independent part of the above two equations to obtain the stationary integral flow 
solution (which is necessary to construct the corresponding phase portrait for the transonic accretion). The time 
independent Euler and the continuity equations are examples of the homogeneous differential equation of first order 
in flow velocity. Integral solutions of such differential equations usually provide the corresponding first integrals of 
motion, - i.e., the quantities which remain conserved. 


For axisymmetric inviscid hydrodynamic accretion we thus have two first integrals of motion obtained by integrating 
the Euler and the continuity equations. 

Since the continuity equation implies the conservation of mass, the corresponding first integral of motion - 
irrespective of the equation of state is used to describe the flow. The integral solution of the Euler equation will 
provide the total specific energy as another first integral of motion for a polytropic flow. 

For isothermal accretion, the corresponding first integral of motion is not to be identified with the specific en¬ 
ergy since energy exchange is allowed for such flow to maintain the space invariance of the bulk temperature of the flow. 


The mass accretion rate, as we will observe, will be a function of the flow thickness, and will resume different values 
for three different geometric configurations of accretion. On the other hand, the first integral of motion obtained from 
the Euler equation does not depend on the flow geometry - it rather depends on the equation of state (though the 

dp 

integration of the term f — contained in the Euler equation). 

P 

b. Adiabatic Flow. — The specific energy first integral of motion can be obtained as 


£ 


u 2 C 2 A 2 

T + Wl + 2e + «’') 


where u is the advective velocity defined on the equatorial plane. 

The mass accretion rates for three different flow geometries are obtained as 


( 12 ) 
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Mch = Hrpu 

(13) 

Mco = Qpur 2 

(14) 


Mve = \ —uCsp 3 ^ 2 (j)' 1//2 

V 7 


(15) 


where H is the half-thickness of the flow and is assumed to be a constant, 0 = H{r)/r , where H(r) is the local 
flow thickness. The corresponding entropy accretion rates can be obtained as 


1 


1 


M = M in j^ - 1 K'y - 1 


(16) 


[71] The subscript [01,00,VE] implies that quantities are evaluated for flow with constant height, in conical 
equilibrium and in vertical equilibrium, respectively. 

The relation between the space gradient of the dynamical velocity and that of the sound speed may be obtained by 
differentiating the expression for the corresponding entropy accretion rates for three different flow geometries. 


dC, 

dr 


= (1-7)' 


CH 


C s (1 du 
2 dr 


u 
2 r 


(17) 


dC s \ C s f 1 du u 

—ZT =( 1 ^T')— o-T- + - 


dr 


co 


2 dr 


(18) 


dCs 

dr 


VE 


1 — 7\ C s 
1 + 7 / u 


du 

dr 


u (3 
2 l r 


(r) 

4>'{r) 


(19) 


We obtain the velocity space gradient (for three different flow geometries) by differentiating the algebraic expression 

dCs 

for £ and by substituting the corresponding values of —— (as obtained by differentiating the respective expressions 


for M) 


dr 


du\ 
dr ) 




CH 


U — 


Cl 


( 20 ) 


CO 




u — 


c? 


( 21 ) 


du\ 
df ) ve 


gi 

7+1 


U — 


2 Cl 
11(7+1) 


( 22 ) 


Following the usual procedure adopted in the literature (see e.g. [14, 58]) the critical point conditions can be 
obtained as 
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Constant Height Model: 


Conical Model: 


Vertical Equilibrium Model: 


u 2 c = C 2 sc = r c (f>'{r c ) - ^ 


n 2 —n 2 = 
yy or U'r 


r c (t>'{r c ) - 


A : 


21 


(23) 


(24) 



<t>'{r c ) + 7 (f>'{r c ) 


( A 2 + r e V(r e ) \ 
\3</>'(r c ) + r c </>” (r c ) / 


(25) 


Substitution of the critical point conditions into the expression for £ will provide an algebraic expression of the 
following generalized form 


£ - /( r cj A, 7 ) = 0 


(26) 


The explicit form of such equations for three different flow models under consideration may be obtained as 
Constant Height Model: 


£- 


1 7 +1 
27 -I 


r c ,CH4>'{r Cl CH) -^— 

r c,CH 


+ 4>(r CtC H) + 


9r 2 

z ' c,CH 


= 0 


(27) 


Conical Model: 


e- 


7 + 1 

7-1 


r c ,co4>'(r c ,co) ~ 


' c,CO 


+ 4>{r c ,co) + 


9r 2 

zr c,CO 


= 0 


(28) 


Vertical Equilibrium Model: 


£- 


[ A2 , vc) ] 

27 

<j>'(r c ,VE) +~f(j>'{r c yE) 

A2 + r l,VE < t>'( r c,VE) 

2 + <P\?cyE) 

Zr c,VE 

7 2 — 1 

r 2 

' C,VE 

3<!>'{Tc,Ve) + ( r c ,VE ) 


= 0 


(29) 


A set of values of [£, A, 7 , a] is required to solve the above algebraic expressions to obtain the value of the corre¬ 
sponding critical point r c . 

£ is scaled by the rest mass energy and includes the rest mass energy itself, hence £ = 1 corresponds to a flow with 
zero thermal energy at infinity, which is obviously not a realistic initial boundary condition to generate the acoustic 
perturbation. Similarly, £ < 1 is also not quite a good choice since such configuration with the negative energy 
accretion state requires a mechanism for dissipative extraction of energy to obtain a positive energy solution [72]. 
Presence of any such dissipative mechanism is not desirable to study the inviscid flow model considered in the present 
work. On the other hand, almost all £ > 1 solutions are theoretically allowed. However, large values of £ represents 
accretion with unrealistically hot flows in astrophysics. In particular, £ > 2 corresponds to extremely large initial 
thermal energy which is not quite commonly observed in accreting black hole candidates. We thus set 1<£ <2. 

A somewhat intuitively obvious range for A for our purpose is 0 < A <4, since A = 0 indicates spherically symmetric 
flow and for A > 4 the sub-Keplerian nature is lost and multi-critical behavior does not show up in general. 

7 = 1 corresponds to isothermal accretion where the acoustic perturbation propagates with position independent 
speed. 7 < 1 is not a realistic choice in accretion astrophysics. 7 > 2 corresponds to the super-dense matter with 
considerably large magnetic field and a direction dependent anisotropic pressure. The presence of a dynamically 
important magnetic field requires the solution of general relativistic magneto hydrodynamics equations which is 
beyond the scope of the present work. Hence a choice for 1 < 7<2 seems to be appropriate. However, preferred bound 
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for realistic black hole accretion is from 7 = 4/3 (ultra-relativistic flow) to 7 = 5/3 (purely non relativistic flow), see, 
e.g., [2] for further detail. Thus we mainly concentrate on 4 / 3 < 7 < 5 / 3 . 

The domain for a lies clearly in between the values of the Kerr parameters corresponding to the maximally rotating 
black hole for the prograde and the retrograde flow. Hence the obvious choice for a is —l<a<l. Although to be 
mentioned here that an upper limit for the Kerr parameter has been set to 0.998 in some works, see, e.g., [60]. We, in 
our work, however, do not consider any such interaction of accreting material with the black hole itself which might 
allow the evolution of the mass and the spin of the hole as was considered in [60] to arrive at the conclusion about 
such upper limit for the black hole spin. 

The allowed domains for the four parameter initial boundary conditions are thus [l<f <2,0 < A<4, , 
4/3<7<5/3, — l<a<l]. 

The four parameter set [£, A, 7 , a] may further be classified intro three different categories, according to the way they 
influence the characteristic properties of the stationary transonic solutions. [£, A, 7 ] characterizes the flow, and not the 
spacetime since the accretion is assumed to be non-self-gravitating. The Kerr parameter a exclusively determines the 
nature of the spacetime and hence can be thought of as some sort of ‘inner boundary condition’ in a qualitative sense 
since the effect of gravity truly requires the full general relativistic framework only out to several gravitational radii, 
beyond which it asymptotically approaches the Newtonian description. [£, A] C [£ , A, 7 ] determines the dynamical 
aspects of the flow, whereas 7 determines the thermodynamic properties. To follow a holistic approach, one needs to 
study the variation of the relevant features of the transonic accretion on all of these four parameters. 

Once the values of the critical points are obtained (physically acceptable values of r c should be validated by the 

du 

criteria r c > r+), transonic accretion solutions can be obtained by integrating the corresponding expressions for — 

dr 

du du 

subjected to the critical value of the — , i.e., the value of — evaluated at the critical point(s). 

dr dr 

Expressions for the critical values of the space gradient of the dynamical velocities can be obtained as (see [18]) 


du 

dr 
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7 (it/) 


(^(r c $'(r c ) - §) + if + $ "(A=)) 


^r c $'{r c ) - 
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(30) 
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(3 4 >"(r c )A 
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.1+7 (&"(r c ) 2 7 /+"'(r c ) 
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c 2 ^$'(r c ) (I+ 7) 2 V $'(rc) 

6 ( 7 - 1 ) <N'(r c ) 6 ( 27 - 1 ) 


7(7 + l ) 2 $'(r c ) 7 2 (7 + l ) 2 


3A 2 


— & (r c ) + — 7 - 


1/2 


(32) 


For certain values of [£, A, 7 , a], one obtains more than one (maximum three) physically acceptable (r c > r + ) values 
of the root of the corresponding 











































£ - /Oc, A, 7) = 0 


(33) 


One can have integral accretion solutions (stationary) with one critical point (saddle type) or three critical points 
where one center type critical point is found in between two saddle type critical points. In subsequent sections, we 
shall develop an eigenvalue scheme to identify which critical point will be of what category. It is to be noted in this 
context that a physical flow can be constructed only through a saddle type critical point. 

At this stage the fundamental distinction between a formal multi-critical configuration and a realizable multi- 
transonic flow requires certain clarification. For stationary axisymmetric accretion, multi-critical flow refers to the 
orientation for which the algebraic solution of the equation expressing the form of the energy first integral of motion 
evaluated at the critical point provides three formal roots for r c all of them being real positive and located outside 
r+. Such configurations, however, does not necessarily produce multi-transonic flows. Critical point behavior is a 
formal property of a differential equation of certain class describing certain behaviors of the first order autonomous 
dynamical systems [61]. On the other hand, transonicity is a physical property which ensures that the flow makes 
a smooth continuous transition for one state of sonicity to the other (from sub / super to the super / sub sonic 
transition). Considering the fact that out of three formal critical points, the middle one being the center type and 
hence not allowing any transonic solutions to be constructed through it, we explain why the realizable multi-transonic 
solutions form a sub-class of the formal multi-critical configurations. 

Subsonic flow starting at the large distance becomes supersonic after crossing the saddle type outer sonic point. Once 
supersonic, flow cannot access another regular sonic point until it is made subsonic again. Certain additional physical 
mechanism (a discontinuous shock transition in the present work). A shock-free multi-critical flow for solution is thus 
a mathematical construct only, for which the in-going accretion solution remains mono-transonic in practice. The 
physically realizable multi-transonic configuration is obtained by allowing the flow solution constructed through the 
outer sonic point with the flow solution passing through the inner sonic point through a steady standing discontinuous 
shock transition. In what follows, we discuss the criteria of such shock formation in somewhat detail. 

c. Shock formation in axisymmetric inviscid accretion. — In this paper we deal with the angular 
momentum supported shock transitions. Certain kind of perturbation may produce discontinuities in large scale 
astropliysical flows. Such discontinuities are supposed to be manifested over a surface when certain dynamical and / 
or thermodynamic flow variables change discontinuously across such surface, and the corresponding surface is marked 
as the surface of discontinuity. Based on specific boundary conditions to be satisfied across the surface of discontinuity, 
such surfaces are classified into various categories - the most important (in astrophysical fluid dynamics) being the 
shock waves or shocks. For adiabatic flow of Newtonian fluid, the boundary conditions satisfied across the shock may 
be expressed as in [62], 


[H] = 0 


(34) 


[\P + Pu 2 }} = 0 (35) 


[[y + fc]]=0 (36) 

where [[/]] denotes the difference of any flow variable / across the shock, i.e. [[/]] = fi ~ fi where [fi, f 2 } are the 
boundary values of the quantity / on the two sides of the shock surface. 

Supersonic astrophysical flows are prone to shock formation phenomena in general and become subsonic after 
the shock. The angular momentum supported centrifugal potential barrier acts as a repulsive agent (against the 
attractive force of the gravity) to break the supersonic incoming flow (matter heading toward a black hole for our 
case). Such shock formation in accretion astrophysics provides an effective mechanism for conversion of a significant 
amount of the ordered gravitational potential energy into ’observable’ radiation by randomizing the inflow. The 
issue of shock formation in this work will be addressed along two different directions. The energy preserving 
Rankine - Hugoniot [62-66] shocks will be considered for the adiabatic flow, where the post shock temperature 
shoots up compared to the position dependent temperature of the pre-shock flow. Higher post-shock temperature 
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puffs up the disc structure which mimics (on local scale) a thick quasi-spherical configuration. The shock thick¬ 
ness is considered to be infinitesimally small - in principle we deal with shocks with shocks with zero effective thickness. 

For isothermal flow we consider temperature preserving shocks. Flow dissipates energy to maintain the invariance 
of the temperature. The post-shock flow thickness does not get altered due to such shocks. 

We define flow variables with the subscript and ’+’ as the ’pre’ and ’post’ shock variables, respectively. For 
constant height flow and flow in conical equilibrium, we obtain the shock condition as 


(7M 2 + l) 2 

2 M 2 + (7 - 1 )M 4 


(37) 


and for flow in vertical equilibrium 


T[M 2 (3 7 -1)+2] 2 T 

_|_2M 2 + (7- 1)M 4 J_ 


(38) 


For isothermal flow, the shock conditions are found to be model independent, i.e., independent of the flow geometry 
and can be expressed as 


p + u+ =p_u_ 


(39) 


P+ + p+u\ = P_ + p-u 2 _ 


(40) 


M * + Ju = M - + TlL 


(41) 
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FIG. 1: Phase portrait of polytropic transonic accretion in conical geometry 


For conical polytropic flow, we show a representative phase portrait for multi-transonic shocked accretion in fig. 1 
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for [£ = 0.0003, A = 3.05, a = 0.5]. ABCDE is the transonic accretion branch through the outer sonic point B. 
The solution JHIJK through the inner sonic point J does not connect the outer boundary of r at infinity. But 
Rankine-Hugoniot standing shock condition connects two branches of solution at CCi and DDi. Hence the accretion 
flow through the outer sonic point may jump discontinuously through DDj to the other branch with higher specific 
entropy being subsonic again and through that branch Di JI it may become again supersonic. The other branch 
FBG in the fig. 1 does not have much relevance in the context of accretion flow. 

The Rankine - Hugoniot condition can be solved to obtain the region of the generic [£, A, 7 , a] which supports the 
formation of steady standing shocks. In fig. 2, we have shown a two - dimensional projection of the four - dimensional 
parameter space (the [£, A] subspace of the [£, A, 7 , a] space for 7 = | and a = 0.5) for conical flow. 



2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 

A 

FIG. 2: The shaded region in £ — A space allowing shock for the conical model 

Region marked by O and I provide monotransonic solutions exclusively passing through the outer type (located 
far away from the black hole) and the inner type (located close to the black hole) sonic points, respectively. Region 
A\ + A 2 represents the multi-critical solutions for which the condition M ln > Ad ou t is satisfied. [£,A, 7 ,cl)a 2 
represents the true multi-transonic region which represents the shock formation. [£, A, 7 , o]a 3 provides multi-critical 
flow with Ad ou t > Adi n . In such region, however, accretion flow remains mono-transonic. 

In fig. 0 c, we plot the [£, A, 7 , a] s hock for three different flow models as shown in the figure. 




FIG. 3: Figure in the left side compares the shock allowing regions for three models. Dashed (blue in online version) lines for 
constant height, dotted (red in online) lines for conical and solid (green in online) lines for vertical equilibrium model. The 
figure in right side magnifies the common region 
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FIG. 4: Variation of shock location with spin for three models at different values of 1 A (see text). Line types and color scheme 
are same as previous figure 


We find that the constant height flow accommodates standing shock for a relatively higher value of A ( A still 
remains sub-Keplerian, though). Fig. Oc shows the region of [£, A, 7 , a] common to all three flow geometries which 
allows the shock formation. 



FIG. 5: Variation of shock strengths with spin for three models (Line types and color scheme same). 


This mutually overlapping shock forming [£, A] region, and its generalized version in the four dimensional [£, A, 7 , a] 
space is very crucial finding of our present work. If one requires comparison between three different models based 
on any physical process / value of the flow variable corresponding to the multi-transonic shocked accretion, one has 
to take [£,A] from the subspace as shown in fig. Oc. In fig. 4, we show the dependence of the shock location on 
the black hole spin parameter for constant height flow (hereafter any variation in the figure corresponding to the 
constant height flow will be drawn using blue color and will be marked by H. For conical flow and flow in hydrostatic 
equilibrium similar variations will be depicted using red and green color, respectively, and will be marked by C and 
V, respectively), conical flow and flow in vertical equilibrium for [£ = 0.0003,7 = 4/3] and for three sets of specific 
angular momentum, A = 2.6 (left most set of curves), A = 3.02 and A = 3.3 (right most set of curves), respectively. 

We find that the shock location non-linearly correlates with black hole spin (same result has been obtained for full 
general relativistic flow in the Kerr metric, see, e.g., Figure 1 in [67]. 
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FIG. 6 : Variation of different shock parameters with a, in comparing three models (line types and color scheme same) 


It is intuitively obvious that the shocks located closer to the black hole event horizon is stronger. This is because 
of the fact that more amount of gravitational potential energy is available to be released for smaller shock location. 
We quantitatively demonstrate such finding through fig. 5, where we have plotted the shock strength (defined by the 
pre to post shock values of the Mach number) as a function of the Kerr parameter. 

Clearly, stronger shocks are produced for smaller shock locations. It is also observed that for the same shock 
location, the shock is strongest for constant height flow and is weakest for flow in vertical equilibrium. In fig. 6 , the 
variation of the ratio of post to pre-shock quantities have been shown as a function of the black hole spin. 

The density compression ratio p+/p-, ratio of temperature T + /T_ and pressure anti-correlate with the Kerr 

parameter for obvious reason (that stronger, hotter and intensely compressed post shock flows are generated closer 
to the event horizon). 

d. Isothermal Flow. — A typical shocked multi-transonic flow topology is shown in fig. 7 for [Tio = 1, A = 
3.05, a = 0.5]. 

The two dimensional projection of the three dimensional [T, A, a] parameter space is shown in fig. 8 for conical 
model with the black hole spin parameter a = 0.5. 

The symbols O and I represent mono-transonic accretion constructed through the outer and the inner sonic points, 
respectively. (A± + A 2 ) represents the multi-critical region for which £,; n > £ out and A 3 represent the multi-critical 
region for < £ out . For Ai and A 3 , however, ingoing accretion flow is mono-transonic, passing through the outer 
sonic point for [T, A] J 4 1 and through the inner sonic point for [T, A] . 4 . 3 . [T, X]a 2 provides the true multi-transonic 
accretion endowed with a temperature preserving standing shock, fig. 9 shows the parameter space spanned by T and 
A for which the shock formation is possible for three different flow geometries. 

Once again, the constant height flow allows a larger [T, A] region for the shock formation compared to other two 
flow geometries. We have a large overlapping region here for the formation of the isothermal shock. The enlarged 
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FIG. 7: Phase plot for isothermal accretion with isothermal shocks, similar to polytropic one (see fig. 1) 



2.6 2.8 3 3.2 3.4 3.6 3.8 4 

A 


FIG. 8: Region of parameter space allowing isothermal shock (see text) 


(sufficiently zoomed in) common region is shown in fig. 9. 

The variation of the shock location with the black hole spin parameter is shown in fig. 10. 

The r s h — a variation follows the same trend as obtained in the polytropic flow. The variation of the shock strength 
with the Kerr parameter has been shown in fig.ll. 

Stronger shocks are formed closer to the horizon and constant height flow produces the strongest possible shock 
variation of the density compression ratio as well as the ratio of the post to pre-shock pressure on the black hole spin 
is shown in fig. 12. 

e. Stability of the stationary flow. — The stability of the stationary flow configuration is to be ensured 
for the validity of this model. Following the method depicted in refs. [14, 58, 68-70], one introduce time dependent 
perturbation of the form u(r,t) = Uo(r) + u(r,t) and p(r,t) = po(r) + p(r,t) where the quantities with suffix 0 
correspond to the values for stationary flow at a particular r and the the additional time dependent parts under 
tilde are the corresponding perturbing parts. These quantities are introduced into the time dependent counterparts 
of continuity equation, eq. (9) and Euler equation, eq. (8). Following the same line of argument as mentioned in the 
ref. [] with general form of pseudo-potential </>, one may show that within the time evolution equation of the perturbing 
part of the mass accretion rate / = purH , the potential does not enter, as long as it being a function of radial 
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FIG. 9: Comparison of regions allowing isothermal shock in parameter space for different models (color scheme same). The 
figure in right magnifies the common region of the figure in left. 



FIG. 10: Variation of shock location with a for three models. 



FIG. 11: Variation of shock strength with a for three models. 
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FIG. 12: Variation of shock parameters with a for three models. 


distance r only. Thus the time evolution equation of the perturbing part / becomes ref. [14] 

^(P^/) =0, (42) 

in which the Greek indices are made to run from 0 to 1, with the identification that 0 stands for t, and 1 stands for 
r; and 



1 u o \ . 
u 0 ul - ac 2 s0 J ’ 


(43) 


a being a model dependent parameter such that a = 1 for CH and CM while a = 2/(7 +1) for VE. Hence following the 
argument of refs. [14, 58], basing upon the identical time evolution equation, it may be concluded that the stationary 
flow is stable in this case too. 

Moreover comparing it with wave equation of a scalar (say ip) in curved space time i.e., 

^7= [fy (V=9<rd v ) !/>] = 0 (44) 

(where g = det(g At[/ )) following [Visserl998], one may form an effective metric for the acoustic perturbation (henceforth 
to be called as acoustic metric, G^ v ) as \J— det(G^)G M1/ = F 1 '. 

/. Acoustic Surface Gravity : Dependence on the black hole spin parameter. — The metric elements 
v corresponding to the sonic geometry is obtained as a consequence of the linear stability analysis. This we have 
demonstrated in the previous section. The acoustic surface gravity n for the stationary background fluid can be 
obtained as (see ref. [18]) 


1 (1 + 2cf>{r)) (1-^-2 <j>(r)^r 


)( 1 

du 

dc s 



/ Vi-** 

dr 

c dr 

C- 

) 


(45) 


As is obvious from the explicit expression of k , the surface gravity is a function of our four parameter initial 
boundary condition governing the flow, i.e., k = k[£. A, 7, a]. 

For a fixed set of values of [£, A, 7], we would like to calculate k at the inner sonic point for a multi-transonic flow for a 
range of Kerr parameters to study the dependence of the acoustic surface gravity on black hole spin, k— a relationship 
could also be demonstrated for the outer acoustic horizon for a shocked flow. However, n at the outer acoustic horizon 
has numerical value way less compared to that of calculated at the inner acoustic horizon. This is a generic property 
(that Ki n >> Kout) found independent of the nature of the background space-time metric, geometric configuration of 
the accretion flow, as well as the thermodynamic equation of state used to describe the matter flow in general. This 
indicates that the numerical value of the acoustic surface gravity correlates with the strength of the gravitational 
attraction of the background gravitational field. Also to mention in this context that n ou t is not sensitive enough on ’a’ 
when evaluated at the outer acoustic horizon. This is intuitively obvious because at the outer acoustic horizon (which 
forms a large distance away from the black hole), space-time becomes asymptotically flat and the effect of the black 
hole spin does not really affect the dynamics of the flow, and hence the nature of the sonic geometry embedded within. 
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FIG. 13: Variation of surface gravity with a for three polytropic flow models. 


In fig. 13 and fig. 14, we show the k — a dependence for adiabatic as well for the isothermal accretion, respectively, 
for three different flow geometries considered in our work. 



FIG. 14: Variation of surface gravity with a for three isothermal flow models. 

Note that the £ — A parameter space responsible for shock formation allows a rather narrow range of the black 
hole spin when we consider the multi-transonic parameter space common to all three flow models. Hence for a fixed 
value of [£ 1 A, 7] or [T, X] , to find an appreciable span of ’a’ is rather non-trivial while comparing the nature of n — a 
variation for different flow configurations. One can, however, bypass such constraints by studying the dependence of 
the acoustic surface gravity as a function of the black hole spin for mono-transonic flow constructed through the inner 
acoustic horizon. For such flow configuration, the entire range of ’a’, i.e., [—1,1], may be accessed and hence a more 
comprehensive information about n — a variation may be obtained. 
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